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oo ' Abstract 
O . 

The nucleoli's axial vector charge, gA, becomes delocalized in the chiral limit. 

O When ttt,^ = 0, and SU(2)lXSU(2)r is exact, 1/3 of the nucleon's axial charge 

Q_i| is to be found at infinite distance from the nucleon. For finite m n this result is 

J|h I approached smoothly as m w — ► 0. We illustrate this effect by considering the 

, lepton-proton spin-spin interaction arising from Z° exchange as a function of 

•i-H . 

^ . m n . Delocalization may have implications for lattice calculations of gA and 

■ in nuclei. 
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In this Letter I point out that the nucleon's axial vector charge, g A , becomes delocalized in 
the chiral limit. Specifically, when m n = 0, exactly 1/3 of the nucleon's axial charge is to be 
found at infinite distance from the nucleon. The remaining 2/3 is found at typical hadronic 
distances. This surprising situation is approached smoothly as m n — > 0: A contribution to 
the axial charge that approaches g A /3 is found at distances of order l/m w when m n is much 
greater than typical hadronic distance scales. |TJ 

To demonstrate that this effect is physical, I compute the effective s-wave spin-spin 
Hamiltonian for a hydrogenic atom composed of a nucleon (mass M) and a charged lepton 
(mass mi) due to Z° exchange, for two cases, a) M ^> m n ^> am f) and b)M> ami 3> m^. 
As expected, in Case b) the interaction is exactly 2/3 as strong as Case a). 

Lattice estimates of the nucleon's axial vector charge extrapolated to the chiral limit 
have always found values below the physical value of g A ~ 1.26. |3]||] These calculations 
are performed relatively far from the chiral limit on relatively small lattices. Although my 
analysis is performed in the unbounded continuum, it can be adapted to a finite domain 
and periodic boundary conditions, and may lead to modification of lattice extrapolations.^] 
Also, estimates of g A in nuclei consistently fall below the value obtained in isolation, [fj] The 
phenomenon discussed here may account for medium dependent modifications of gA without 
invoking partial restoration of chiral symmetry. 

At first I work in the chiral limit, where the pion is massless and SU(2)l x SU(2)r is an 
exact symmetry. After deriving the claimed result, I present a simpler, heuristic argument 
based on symmetries, which may help develop the generalization to other geometries. Next, 
I consider the situation when the pion is light but not massless, and show that the chiral 
limit is approached smoothly. Finally I mention extensions to other hadronic axial charges 
and point out that this anomaly is special to the axial current and does not afflict either 
the nucleon's tensor charge or the matrix elements of the tower of twist-two operators that 
determine its polarized quark distribution, Aq(x,Q 2 ). 

The nucleon's axial charge is defined by the matrix element of the isovector axial current 
between momentum eigenstates: || 



(N(p', s')\A(0)\N(p, s)) = u(p', s')\t [g A (q 2 )ll, + h A (q 2 )q l5 ) u(p, s). (1) 

The axial charge, g A , is <7a(0). h A (q 2 ) is the induced pseudoscalar form factor, which receives 
a contribution from the pion pole, 

7 / 2 \ d A (q 2 ) 2f n g nN 

hA{q = — = — + ■ ■ • 2 



1 For an analysis of the implications of this paper for lattice calculations see Ref. B. 
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where g w N is the 7r-nucleon coupling constant, f w is the pion decay constant, and the terms 
denoted by • • • are less singular than 1/q 2 as q — > 0. Since the isovector axial current is 
conserved in the chiral limit, g& and g^N are related by the Goldberger-Treiman relation J7| 

9A = -j^- (3) 

which is exact when m n = and approximate for small 

Consider now the matrix element of the axial current in a nucleon wavepacket state, 
|iV(s)). N(s) represents a nucleon with spin s at rest, centered at the origin. It is de- 
scribed by a wave packet normalized to unity, \N(s)) = J , d J/ 2 0(p)|-^(p, s)), (N(s)\N(s)} = 
J d 3 pjj\(f)(p)\ 2 = 1. Consider the fourier transform of the current matrix element in the limit 
where \q | is small compared to the momentum scale of the wave packet. Using eq. ([!]) and 
the following properties of Dirac spinors, 

u(p + q, s)-f 5 u{p, s) -> T^V^iPi s h»l5u(P, s) as g M -> (4) 



I obtain 



lim 



ay tfVe^(iV( S ) |4(f) |AT( S )) = ^(^-g ( 5 ) 

up to recoil corrections of order (p 2 )/m 2 . \ encodes the spin-isospin content of the matrix 
element, 

X = \ularU Q . (6) 

Uq is a spinor-isospinor for a nucleon at rest. In deriving eq. (||) I have used the Goldberger- 
Treiman relation to express the entire matrix element in terms of g&. Conservation of 
the axial current in the q — > limit is manifested by the fact that q contracted with 
eq. dB|) vanishes. In this calculation it is necessary to treat the q — > limit carefully. Since 
q° rv> q 2 /2M for small q, I have replaced q 2 by — q 2 . Also from eq. (f|), u(q, s)7 5 u(0, s) — > 
2J/?7qC? • q Uq as <f — > 0. The term proportional to q q ■ a in eq. (^) shows that the limit 
g — > is nonuniform. 

To separate the short- and long-range contributions to the matrix element I define con- 
tributions to eq. (§) from within and outside a sphere of radius R centered on the nucleon. 
R is taken much greater than the range of typical QCD interactions, R^> 1/A, and is held 
fixed as q — > 0, 



2 Other form factors, which are irrelevant for this discussion, have been omitted from eq. (|l|). 
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(N\A(r)\N) = 6{R - r) a (r) + 6{r - R) a n (r). (7) 

The first term is the unspecified contribution to the axial current from typical hadronic 
distance scales. For large enough R only the single pion contribution to the axial current 
remains when r > R, so the second term is completely due to the pion, A — ► — /^V^, which 
in turn is given by the familiar pole diagram. 

Consider first the contribution of 8(R — r) a^. Because the integral is bounded by r < R, the 

q — > limit can be taken without any special precaution. Therefore this piece of the axial 
current can only contribute to g A , not to the pion pole. Thus this contribution to eq. (||) is 
given by 

i-R , .. ... /■ <" 

9^0 J 



lim [ R r 2 drdtt e iqr %(r) = r 2 drdtt a {r) = g° A {R)x- (9) 
q^oJo Jo ~ 



There is no need to specify g\{R) further, except to emphasize that it is the result one would 
obtain by restricting the calculation to the interior of a sphere of radius R. 

Next turn to the region r > R, where the pion contribution to the axial current is explicit 
and exact, 

J™ r 2 dr j dVL e^%(r) = -C J™ r 2 dr j dVt e** V(— r ■ x) (10) 

where C = g^j^jAn^M. It remains merely to evaluate this integral, mindful that q — > 
with R fixed. 

Note that dimensional analysis suggests that the contribution from eq. fllOD will persist 
as q — > no matter how large R is, because the r integration is scale invariant. It must 
generate the pion pole term, — q q-X9A in eq. (|5[). The surprise is that there is another term 

proportional to x alone, which contributes to g A - 

Eq. (|It]) may be evaluated directly at fixed R and q. In particular, the oscillatory q 
dependence ensures convergence at large r and does not require any special treatment. It is 
somewhat easier however to integrate by parts: 

R 

rco r roc r r 

/ r 2 dr & dtt e^%(f) = Ciq / dr & dtte iq ^ r-x+C & dtte iq ^r f-\ 
Jr J Jr J ~ J ~ 

= r v + r R + r 00 , (11) 

where r y is the volume contribution for r > R, and T # i00 are surface terms at R and at 
the surface at infinity. 
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The surface at infinity does not contribute. To see this, expand the exponential using 
exp(iq ■ r) = ^2ii i (2£ + 1)-Fje(g • f)ji(qr). From the physical definition of g& it is clear that 
T oo is to be evaluated on the surface at infinity with q and R fixed. Only I < 2 contributes 

and the Bessel functions ji(qr) go to zero as r — > oo at fixed q, so Too = 0. 

Next consider the volume contribution, f ■ x selects out the £ = 1 term in the Bessel 
expansion of the exponential: 

— > Z" 00 f 

Tv = —3Cq J drji(qr) J dVl q-f f-\ 

= -AnCq q-x jo(qR) (12) 

where I have used {r^j) = Taking q — > at fixed i?, I obtain 

i- -p 9-rrNfir » * ^, /io\ 

hml y = — — o o-y 13) 

which is the expected pion pole term. 

Finally consider the surface term at R. With R fixed and q — > 0, the exponential can be 
replaced by unity with the result 



limT b = C dVt f r-v 



3M 



X (14) 



which contributes to qa- Using the Goldberger-Treiman relation the magnitude of this term 
is qa /3 independent of R. 

The three contributions to eq. fl5|) given in eqs. @, (0), and ( jTjj) combine to give 



hm / d 3 re^(N(s)\A^\N(s)) = g A (R)x+\9AX-9A q q~X (15) 

again using the Goldberger-Treiman relation. 

Comparison with eq. fl5|) verifies that we have obtained the correct pion pole contribution 
(a check on the algebra) and that the contribution to Qa from inside an arbitrary but fixed 
radius R equals |(yu independent of R. 

Since this result is unexpected it is worth presenting a somewhat different argument, 
which is simpler, but heuristic. This argument may generalize to other geometries. If one 
suspects a possible long range contribution to the matrix element on the left-hand side of 
eq. (H) one can suppress it by averaging over the angles of q before sending q — > 0. This has 
the effect of replacing exp(iq-f) by sin qr/qr, which is suppressed by an additional power of r 
at large r. Thus the spherical average of eq. (|5]) is equivalent to the short-range contribution 
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piece of eq. fl5|) and should be identified with g\(R) X- A factor of 1/3 appears when one 
averages q q-a over the sphere, with the result 

Hm<y d 3 re^(N(s)\A(r)\N(s))^ = g° A (R)x=l9AX (16) 

which reproduces our earlier result. 

Thus I conclude that in the chiral limit the correct value of the nucleon's axial charge is 
3/2 the value that is computed by restricting the calculation to any finite volume. 

Next, I repeat the analysis with a nonzero pion mass and study the behavior of g A as a 
function of m n and R. Although the form of the analysis changes considerably, the physical 
result changes smoothly as one departs from the chiral limit. There is no longer a pion 
pole in the induced pseudoscalar form factor, so that term plays no role in the analysis as 
q — > 0. Note also that the Goldberger-Treiman relation is no longer exact so I will not 
replace g W Nf n /M by g A except in the limit m w — > 0. All contributions to the matrix element 
of A(r) fall exponentially, ~ e~ m * r , for large r. Still, for small m v and large r (m^, l/r< A) 

the pion dominates and its contribution can be calculated analytically. In place of eq. (H) 
one has 

9nNfn, 



a w (r) 



A-kM 



VrV(V)' ( 17 ) 



It is no longer necessary to be so careful about the limit q — > 0. Instead the matrix element 
of the axial current at q = can be calculated directly. It can be separated into the 
contribution from r < R, again labelled g A (R), and the contribution from r > R where the 
pion dominates. A straightforward calculation yields 

d 3 r(N(s)\A(r-)\N(s)) =9aX= (g° A (R) + l^f^i 1 + ™^)e" m ^) X ■ (18) 

This result displays the expected behavior: For any fixed R, as — > 0, one third of the 
axial charge comes from outside R. However, for any fixed m n , as R — > oo all the axial 
charge comes from inside R. Finally, it is possible to calculate the density of axial charge in 
the region where the pion dominates. The result, 

dgA _ l fhNf* 2 m , r , v 

dr ~ 3 M * 6 1 9J 

valid at large r, again displays the nonuniformity of the —>■ and r — > oo limits. Whether 
these results are modified by chiral logarithms for small but nonvanishing is outside the 
scope of this Letter. 

A comment is in order on the relation of this result to traditional calculations of the 
pionic contribution to the nucleon's axial charge. M It is well known that, when integrated 
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over all space, the pion contributes to /?-a(<? 2 ) and not to gA(q 2 )- This can be seen by Fourier 
transforming eq. (|T7| ) with the result, 



(assuming some smearing of the pion source at the scale of the nucleon size). The tensor 
structure qX'Q signals a contribution to liA(—q 2 ) and no contribution to <m(— <? 2 ). There 

is no puzzle here: we are studying the contribution to the axial charge outside a sphere of 
radius R surrounding a localized nucleon, not the integral over all space. The result we 
have derived is independent of the nature of the short range contributions to the nucleon's 
axial charge and the fact that the total pion pole contribution vanishes. Finally, the reader 
may wonder what additional short range term appears in the pion contribution to dgA/dr 
(see eq. (0)) to make the integral over r vanish. A study of eq. (O) at small r reveals the 
presence of a Dirac 5 function, so the pionic contribution to dgA/dr is 



which integrates to zero. For an extended nucleon source, presumably the 5 function would 
be replaced by a smooth contribution with the same integral. 

One might worry that the derealization of the axial charge demonstrated here is a merely 
a mathematical oddity rather than a physically observable effect. To put this doubt to rest, 
consider the following gedanken experiment in which it could, in principle, be measured. 
Consider a hydrogenic atom composed of a proton bound to a lepton of mass rri£. We 
take me small enough that the Bohr radius, ai = l/m^a, is much greater than 1/Aqcd, 
the range of the strong interactions, excluding the pion which we treat explicitly. Consider 
the proton and lepton interaction mediated by Z° exchange. The axial-axial piece of Z° 
exchange generates a spin-spin interaction which contributes a very small (but in principle 
observable) correction to the hydrogen hyperfme splitting. An elementary calculation yields, 




(20) 




pion 



(21) 




G p 



(22) 



where \P(s)) is a proton with spin s localized at the origin, A% is the third component of A 

in isospin space, and \<$i is the lepton spin operator. We are interested in the chiral limit, 
so we take m n <C A QCD < M. We parameterize A in terms of a and a n (see eqs. (^) and 

(|T7|). In the s-wave only the spherical average of A3 survives, and we find, 




(23) 
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— * — * 

where ((A 3 (r))) denotes the spherical average of the proton matrix element of A 3 (f). I 
have replaced the short range contribution, ao(r)8(R — r) by g\x^{^) since the hydrogenic 

ground state wave function is essentially constant for r < R. Also, in averaging the pionic 
contribution I have omitted the delta function at the origin in V 2 (e _m7rr /r) because the 
origin is explicitly excluded by the 8(r — R). Combining eqs. fl22| ) and (^), yields an 
effective Hamiltonian for the spin-spin interaction in the s-wave, 

G { 1 1 e — m " !rr \ 

#s-wave(r) = l-g A 5 3 (r) + ^9Avn\—^— \ o P ■ a e (24) 



Next, take the expectation value of if s _ wave (f) in the ground state ipo(r) = (^) 3//2 e r ^ ai . If 
l/m-x ^> ae, then the second term in eq. ( p3|) is negligible and we obtain, 



G fl \ 

lim (^ol-^s-wavel^o) = ( ~ 9a J &p • <^l^o(0)| 2 (25) 

which verifies that one-third of the nucleon's axial charge has receded to distances of order 
l/m w , which is much greater than the size of the atom and therefore does not affect the 
interaction. On the other hand, if l/m n <C then the second term in eq. fl23|) is easily 
evaluated and supplies the missing third of g A - 

G 

lim (V'ol-ffs-wavelV'o) = -—f=g A o P ■ a e \i> (0)\ 2 (26) 



m-n-a^— > oo 



The latter case goes smoothly over to the physically relevant case (with strength g A ) in 
which m n and Aqcd are comparable and much greater than This example shows that 

the derealization of the axial charge in the chiral limit is (at least in principle) a physically 
measurable effect. 

Other hadron axial charges, for example, the ones that enter strangeness conserv- 
ing semileptonic hyperon decays (like £ — > Aei>) or semileptonic D and B decays (like 
B* — > Beu), behave in the same way as m n — > 0. However, it appears that other nu- 
cleon matrix elements of pseudotensor operators do not suffer from this subtlety that afflicts 
the axial charge. The isovector tensor charge, gx, defined as the coefficient of p fJj s u — p v s^ 
in the forward matrix element of ipo'^'j^^Tip, provides an important contrast because it 



3 Note that in this calculation it is justified to ignore the q 2 variation of the axial form factor 
because the typical momentum transfer in the hydrogenic ground state is determined by m^, which 
can be taken to be small compared to the natural scale of variation of the axial form factor (which 
is finite in the chiral limit). |§] Note also that I have identified g-KNfn/M with g A even for 7^ 0, 
which is only approximately true. 
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is amenable to lattice simulation. There is no pion pole term in the tensor charge be- 
cause (7r(p) 1^(7^75 1 Ttp\0) = 0. So lattice calculations of gr should be more stable in 
the chiral limit than qa- Another important tower of operators are the twist-two, pseu- 
dotensor operators that determine the spin-dependent quark distributions in the nucleon, 
•A%,iH2-nn = V'7{mi-^M2 ' ' ' D^^ip. A n does receive a pion pole contribution, but it is pro- 
portional to a high (second or greater) derivative of 7r, and therefore does not generate a 
nonuniformity in the matrix element. 

It appears that the flavor nonsinglet axial charge is unique in possessing such singular 
behavior in the chiral limit. This may be the source of instabilities in lattice simulations of 
qa- The unusual spatial distribution of the nucleon's axial charge may have other physical 
consequences. For example, that it receives an important contribution from distance scales 
of order l/m n may account for modifications of qa in large nuclei without invoking partial 
restoration of chiral symmetry. 
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